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Abstract. We study a method to determine the residual conductance of a correlated system by means of 
the ground-state properties of a large ring composed of the system itself and a long non-interacting lead. 
Lhe transmission probability through the interacting region and thus its residual conductance is deduced 
from the persistent current induced by a flux threading the ring. Density Matrix Renormalization Group 
techniques are employed to obtain numerical results for one- dimensional systems of interacting spinless 
fermions. As the flux dependence of the persistent current for such a system demonstrates, the interacting 
system coupled to an infinite non-interacting lead behaves as a non-interacting scatterer, but with an 
interaction dependent elastic transmission coefficient. Lhe scaling to large lead sizes is discussed in detail 
as it constitutes a crucial step in determining the conductance. Furthermore, the method, which so far had 
been used at half filling, is extended to arbitrary filling and also applied to disordered interacting systems, 
where it is found that repulsive interaction can favor transport. 

PACS. 73.23.-b Electronic transport in mesoscopic systems - 71.10.-w Lheories and models of many- 
electron systems - 05.60.Gg Quantum transport - 73.63.Nm Quantum wires 



1 Introduction 

Large experimental activities have recently been de- 
voted to the study of the conductance of low-dimensional 
nanosystems like molecules, atomic chains, nanotubes, 
and quantum wires [T||2||^||4"||^] with sizes typically of the 
order of the electronic Fermi wavelength. Since the screen- 
ing of the Coulomb interaction in such systems is less effec- 
tive than in three dimensions, electronic correlations can 
no longer be neglected with respect to kinetic effects. In 
some of the systems mentioned, the Luttinger liquid be- 
havior [S][7| is relevant and might influence the transport 
properties. 

The correlations become particularly relevant for low 
temperature electronic transport properties like the resid- 
ual conductance and the interpretation of the experi- 
mental data requires a good understanding of transport 
through a region with strong correlations. However, this 
turns out to be a demanding task and various attempts 
have been made in this direction 001]. 

The purpose of the present work is to contribute to 
the fundamental problem of transport through correlated 
nanostructures by studying a novel approach where the 



conductance is obtained from thermodynamic properties 
of a ring consisting of the nanosystem and a long lead. 
Such an embedding method has been actively pursued in 
the last few years |1()U1 111 12111. 'ill 1 111 15111 (ij . Here, we crit- 
ically study its hypotheses and consequences in order to 
put it on a firm theoretical basis. 

A powerful concept which was used for studying co- 
herent transport through non-interacting systems is the 
Landauer-Biittiker formalism |18IIL9| which formulates a 
scattering problem between electron reservoirs. Although 
the electrons in the reservoirs interact, their density is 
very high such that the Coulomb energy to kinetic en- 
ergy is small and they can be replaced by non-interacting 
quasiparticles. Hence, the reservoirs are well described by 
a Fermi distribution characterized by a temperature and 
a chemical potential. Within the scattering approach, the 
dimensionless residual conductance g (in units of e 2 /h) is 
given by the elastic transmission probability \t(E-p)\ 2 at 
the Fermi energy Ep . 

The situation becomes more complicated if electron- 
electron interaction is present in the scattering region be- 
cause the passage of electrons may lead to the creation of 
excitations. However, for temperatures smaller than the 
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Fig. 1. The system considered within the embedding approach 
is a one-dimensional ring consisting of an interacting region 
(grey) of length Ls and a non- interacting lead (black) of length 
Ll. The ring is threaded by an Aharonov-Bohm flux cf>. 

characteristic excitation energy of the nanosystem, the 
idea of the Landauer-Biittiker formalism still applies |S] 
because all accessible states in the reservoir with an en- 
ergy lower than the excitation energy are occupied. In- 
elastic processes are then forbidden. On the other hand, 
it remains non-trivial to determine the clastic transmis- 
sion probability through a correlated system. Green func- 
tion methods, while being conceptually adequate, require 
knowledge of the excited states and may become numeri- 
cally quite involved. 

An alternative approach consists in considering the 
ground state properties of a ring formed by the system 
of interest, which we will refer to as correlated system or 
nanosystem, and a very long non-interacting lead as de- 
picted in Fig. ^ Within this embedding method, the rele- 
vant information about the conductance can be extracted 
by means of a flux threading the ring, which gives rise 
to a flux dependence of the ground-state energy and thus 
to a persistent current. This setup accounts for two im- 
portant physical ingredients of coherent transport. First, 
the flux dependence of the ground-state energy provides 
information about extended states in the interacting re- 
gion. Second, the two contacts between system and lead 
allow to transfer electrons into the system. This is an es- 
sential point in the description of conductance |2U| , which 
is not present when the persistent current is calculated for 
a correlated system without auxiliary lead. 

Favand and Mila used the above described approach to 
compare, within a model of spinless fermions, the tunnel- 
ing conductance of molecules with a Mott-Hubbard gap 
and of molecules with a dimerization gap |1U) . Sushkov 
used the same idea for a study of the 0.7e 2 /h anomaly 
observed in quantum point contacts |1 1U21| . However, an 
important difference with respect to Ref. ^U] is that he 
kept the interaction in the leads within the Hartree-Fock 
approximation. As the present authors have emphasized 
|12| . the extrapolation to infinite lead length can only yield 
meaningful results if no interaction is present in the aux- 
iliary lead. Other important aspects discussed in Ref. ^2 
are the relevance of the contacts, the oscillation of the 
conductance as a function of the number of sites in the 
interacting region, and the role of static disorder. Meden 
and Schollwock compared the results obtained within this 
approach to those of a perturbative functional renormal- 



ization group and showed that both give the same results 
at small values of the interaction strength, verifying scal- 
ing laws associated with Luttinger liquid behavior |13U14j . 
Rejec and Ramsak tested the method, comparing its pre- 
diction with previous results for transport through single 
and double quantum dots. They presented a generaliza- 
tion to systems without time-reversal symmetry, using as 
an example a nanosystem which itself forms an Aharonov- 
Bohm ring |15U16j . 

An approach related to the embedding method has re- 
cently been proposed by Chiappe and Verges ^2] m which 
the nanosystem and a small part of the leads are diagonal- 
ized exactly. In a second step, this subsystem is attached 
to semi-infinite leads and Green functions arc employed to 
numerically calculate the conductance. The conductance 
through a one-dimensional interacting spin-system cou- 
pled to non-interacting leads was also studied by Louis 
and Gros by means of a Monte-Carlo based method \'22\. 

The relationship between the conductance and the per- 
sistent current of a large ring has only been proven for 
non-interacting scattering systems. No rigorous proof has 
so far been put forward once electronic correlations are 
present in the scattering region. However, the conductance 
obtained by means of the embedding method satisfies all 
basic requirements and reproduces the correct behavior in 
various limiting cases. Moreover, in this work we demon- 
strate numerically for the one-dimensional case that in 
the limit of very large ring size, the effect of an interact- 
ing scatterer on the persistent current can be described 
by the amplitude of a transmission probability character- 
izing a 2 x 2 transfer matrix. Thus, transport through an 
interacting region can be understood as a non-interacting 
scatterer with interaction dependent parameters. 

The remainder of the paper is organized as follows. In 
Section [21 we use Density Matrix Renormalization Group 
(DMRG) techniques to calculate the flux dependence of 
the persistent current through a ring composed of an in- 
teracting region and a non-interacting auxiliary lead in 
the limit where the latter becomes very long. It is found 
that this flux dependence reproduces the one expected for 
a non-interacting ring of equal length interrupted by a 
scatterer which can be characterized by a transfer matrix. 

In the absence of Luttinger-like correlations in the ring 
it is meaningful to consider the limit of a very long auxil- 
iary lead. In Section |3| we will explain how the extrapola- 
tion to infinite circumference can be performed in order to 
extract the interaction-dependent transmission coefficient 
and thus the conductance. For this scaling analysis, we 
make use of the charge stiffness instead of the persistent 
current, because it provides us with the same information 
but requires less numerical effort. Specific attention will 
be paid to the case of resonances, which appear when the 
coupling between system and leads is small and where the 
extrapolation has to be done with particular care. 

In the literature, the embedding method has so far 
been discussed only for the case of half filling. In Sec- 
tion 01 we will present an extension to arbitrary filling. 
The important point is to choose the appropriate com- 
pensating background potential which ensures the correct 
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charge density in the system even in the presence of in- 
teractions. While at half filling, it is straightforward to 
define the compensating potential from particle-hole sym- 
metry, a self-consistent procedure is required away from 
half filling. In Section [SJ we employ this new method to 
demonstrate that strong repulsive interactions can favor 
zero-temperature transport through strongly disordered 
systems. 

We present our conclusions and perspectives in Sec- 
tion In Appendix ^ we address the flux dependence 
of the ground state for a ring containing a local non- 
interacting scatterer, and obtain the asymptotic values 
and the finite size corrections to the charge stiffness. In 
Appendix IbI we extend the approach to superconducting 
nanosystems and verify that the known behavior resulting 
from Andreev scattering at the two extremities of a super- 
conducting nanosystem is reproduced. This illustrates the 
validity of the studied embedding method in an extreme 
limit where an attractive electron-electron interaction has 
dramatic effects. 



2 Flux dependence of the persistent current 
for large rings with a small scattering region 

The aim of this section is to demonstrate that the trans- 
port properties of an interacting region can be described 
as a non-interacting scattering problem with interaction 
dependent parameters. We start by considering the setup 
shown in Fig.^which will be employed to study the trans- 
port properties of a one-dimensional system of length L§ . 
This system may contain a scattering potential and, pos- 
sibly, electron-electron interaction may be present there. 

The system is contacted by the two ends of an aux- 
iliary one-dimensional lead of length Ll so that a ring 
of total length L = L$ + is formed. From this setup, 
transmission properties of the system can only be deduced 
if Luttinger liquid correlations |2.'{j in the one-dimensional 
ring are absent. Therefore it is crucial that in the auxil- 
iary lead no electron-electron interaction may be present. 
Not only, this allows to avoid Luttinger liquid correla- 
tions, but the electrons of the combined ring form a Fermi 
liquid in the limit of infinite lead length. According to 
Sushkov, one can give a general argument for Id spinless 
fcrmions on a ring demonstrating that they form a Fermi 
liquid though interactions act in a region of the ring , 
as far as it remains finite while the non-interacting lead 
becomes infinite. This is corroborated by our numerical 
findings presented below. 

Information about the transmission amplitude |£(-Ef)| 
at the Fermi energy E F can be obtained by means of 
a magnetic flux <fi threading the ring. For convenience, 
we introduce the dimensionless flux <P = 2ir(f>/<fio where 
tpo = h/e is the flux quantum. The many-body ground 
state energy Eq of the ring will oscillate with period </>o 
as a function of the flux. The magnetic flux threading the 
ring breaks the symmetry between left and right moving 
electrons and thus gives rise to a persistent current J(<£), 
which at zero temperature is given by J(^) = —dEo/d(f>. 



For non-interacting scatterers, the persistent current 
J(<?) decreases like 1/L for large circumference L of the 
ring. The leading contribution is found to read |25| 



_ ev F Arccos(|i(fc F )| cos(#)) 



Arccos(|t(fc F )|cos(^)j )]sinm (1) 
irL ^1- |i(fc F )| 2 cos 2 (<Z>) 
for an odd number of particles and 



ev F Arccos( \t{k F ) cos £ - tt ) . 
— j 71 ;== = == \t{k F ) sui($) 2 

tt£ yj\ - \t(k F )\ 2 cos 2 (<£) 

for the case of an even number of particles in the ring. 
By Arccos, we denote the principal branch of the inverse 
cosine function which takes values in the interval [0,7r]. 
The derivation of these results is outlined in Appendix lAl 
The persistent currents Q and depend on the 
properties of the non-interacting scatterer only through 
its transmission probability \t(E F )\ 2 at the Fermi energy. 
This important feature allows us to determine the trans- 
mission probability and thus the residual conductance 
of the system from the persistent current of the com- 
posed ring. The relation becomes particularly simple for 
<P = 7r/2, where the transmission coefficient at the Fermi 
energy can be expressed as |2l)lll()llll| 



|t(£ F )f 



W (tt/2) 



(3) 



Here, J° is the persistent current for a clean ring of length 
L. 

We now turn to an interacting nanosystem and demon- 
strate numerically that, in the limit of an infinitely long 
lead, the flux dependence of the persistent current is of the 
same form as in the non- interacting case of Eqs. Q and 
(J2J. The interaction thus enters the result only through 
the transmission coefficient \t(E F , U)\ . 

Specifically, we have performed direct numerical calcu- 
lations of the persistent current for a tight-binding model 
with N interacting spinless fermions on L sites described 
by the Hamiltonian 

L L s 

H = -*E(4<*-i+4-iO+y>[n* " V+] [nt-i - V+] . 



i=l 



i=2 



The hopping amplitude t between nearest neighbors will 
be set to 1 and thus defines our energy scale. c i (cj) is the 
annihilation (creation) operator at site i, rii = c\c i is the 
number operator, and the flux enters through the bound- 
ary condition cq = cxp(i<P)cL. The length scale is given by 
the lattice spacing and the interaction acts between near- 
est neighbors inside the sample (sites i = 1 to Ls), but 
vanishes in the lead. To avoid depletion of electrons in the 
sample due to the repulsive interaction, we introduce a 
compensating potential V + that acts as a positive back- 
ground charge and ensures the local charge neutrality. For 
a half-filled ring, the compensating potential is equal to 
the filling factor v — N/L. Thus, 



v = - = - 
2 J 2 



(5) 
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Fig. 2. The scaling of the persistent current with the total 
length L is performed for a ring with system size Ls = 6 at 
half filling for several values of the flux and the interaction 
strength U. The persistent current J(^) is depicted for even 
particle numbers TV = 6, 8, 10 while for odd particle numbers 
N = 7,9, results for J(<£) with $ = <P — n are shown. The 
extrapolation L — > oo has been performed by means of fits to 
second-order polynomials in 1/L. 

guarantees particle-hole symmetry even in the presence of 
interactions. Outside half filling, this symmetry is broken 
and the compensating potential V+ becomes a function of 
U , N , Ls and Ll as we will discuss in Section 

For the model J2J with system size Lg = Q and at half 
filling, we have numerically determined the persistent cur- 
rent </(<£) supported by the ground state for various values 
of Ll by means of a complex DMRG algorithm. With this 
implementation, we are able to treat not only the flux 
values <P = and <P = ir used in where the Hamilto- 
nian Q can be represented by a real matrix, but also the 
general case of arbitrary flux where the matrix becomes 
complex. In order to determine the persistent current, we 
directly evaluate the current operator for the ground state, 
thereby avoiding the potentially difficult procedure of tak- 
ing numerically the derivative of Eq(<P). 

The length dependence of the persistent current and 
the extrapolation to infinite lead length is shown in Fig. [21 
for particle numbers N = L/2 between 6 and 10. Moti- 
vated by the symmetry 

J(<Z>; N odd) = J(<2> - tt; N even) , (6) 

valid in the non-interacting case according to QJ and (J5J) , 
we plot the interacting results corresponding to even and 
odd N at flux values <P and <P — tt, respectively. As is shown 
in Appendix A, the scaling laws for even and odd N may 
be different. However, making only the flux transformation 
of Eq. allows us to obtain good asymptotic results from 
a single fit to the ensemble of data points for even and odd 
N. A second-order polynomial fit describes very well the 
deviation of the logarithm of the persistent current from 
its asymptotic value. 



Fig. 3. The flux dependence of the persistent current for a 
system size Ls = 6 and half filling is shown for interaction 
strengths U — 1 and 4. The points represent DMRG results 
extrapolated to the limit of infinite leads (see Fig. [5Jl. The 
lines represent the theoretical result for a ring with a non- 
interacting scatterer and transmission amplitudes \t\ = 0.938 
(solid line) and |t| = 0.425 (dotted line). 



The results presented in Fig. [21 indicate that the sym- 
metry (JHl holds even in the presence of electron-electron 
interaction and is independent of the interaction strength 
U. This provides numerical evidence that it should be pos- 
sible to relate the persistent current in the presence of 
an interacting region to the persistent current for a non- 
interacting scattering problem. 

The flux dependence of the persistent current </(<£) for 
an even number of particles, extrapolated to the limit of 
an infinite lead, is presented in Fig. [31 for moderate and 
strong interaction, U = 1 and U = 4, respectively. At the 
filling factor v = 1/2 used here, the interaction effects are 
expected to be most important. As can be seen from Fig. EH 
the flux dependence of the persistent current is described 
very well by the expression ^ for the non-interacting case 
with transmission amplitudes of \t\ — 0.938 (solid line) and 
0.425 (dotted line) for U = 1 and 4, respectively. 

This demonstrates that, in the limit L — > oo, the zero- 
temperature persistent current of a ring containing an in- 
teracting region is quantitatively described by the persis- 
tent current of a ring with a scatterering region. A single 
parameter, the interaction-dependent elastic transmission 
coefficient at the Fermi energy \t(Ep, U)\ 2 suffices to char- 
acterize the interacting sample, at least as far as the flux 
dependence of the ground state energy at zero tempera- 
ture is concerned. 

We emphasize that the DMRG technique employed 
here to calculate the persistent current of the ground state 
of the Hamiltonian Q does not rely on any assumption. 
In particular, the DMRG technique does not require that 
the correlated nanosystem must be a Fermi liquid. But 
the fact that the expressions (JIJ and Q for the persistent 
current hold in the infinite lead length limit even in the 
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presence of an interacting region provides strong evidence 
that the Fermi liquid behavior is retained in this limit. 
This result is in agreement with the theoretical expecta- 
tion mentioned above. 

Our findings constitute a numerical "proof" that the 
extension of the relation between persistent current and 
transmission from a non-interacting to an interacting sys- 
tem is correct. Assuming that the composed ring forms 
a Fermi liquid, a discussion of the relation between the 
persistent current and the conductance had already been 
given in |16| . Together with the results of this section, this 
opens a road towards the calculation of the conductance 
for interacting nanosystems. 



3 Conductance from transmission for 
interacting scatterers 

Instead of the persistent current, we will, in the following, 
mostly work with the charge stiffness defined as: 




L 



D = (-1) N -{E(0)-E(n)) 



(7) 



Fig. 4. The scaling of the logarithm of the charge stiff- 
ness with the ring size is shown for systems at half filling 
and Ls = 20, U = 3 (squares), L$ = 12, U = 4 (diamonds), 
Ls = 17, U = 1 (circles), and Ls = 13, U = 2 (triangles). 
The lines are linear fits to the large- L behavior, providing the 
extrapolation to infinite ring size, i.e. 1/L — > 0. 



which describes the change of the ground-state energy 
from periodic to antiperiodic boundary conditions. The 
factor (— 1)^ renders D positive because the many-body 
ground state is diamagnetic for odd N while it is param- 
agnetic for even N. This fact was proven by Leggett 
for spinless fermions in the presence of arbitrary one-body 
potentials and arbitrary strength of electron-electron in- 
teractions. We prefer to work with the charge stiffness D 
instead of the persistent current J because it allows to 
avoid the use of a complex implementation of the DMRG 
algorithm and thus reduces the numerical effort. 

For the case of a non-interacting scatterer, the flux 
dependence of the ground-state energy is derived in Ap- 
pendixE] ^From Eqs. and J2EJ it follows that for the 
limit of infinite lead length we have 



D = 



Arccos(|£(fcF)|) 



(8) 



independent of the parity of N. Solving (JHJ for the tun- 
neling probability yields f^2] 



\t{k F )\ =sm[ - — 



(9) 



where D° is the charge stiffness for a clean ring of length 
L in the absence of electron-electron interactions. We note 
that for weak transmission (\t\ <C 1), D is proportional to 
\t\. 

We have verified that the transmission coefficients cal- 
culated from the stiffness using Eq. as described in 
Ref. ^2] coincide with the ones obtained by fitting the 
full flux dependence of the persistent current (Fig.|3J to a 
precision better than 0.5%. 



3.1 Scaling of the stiffness and extrapolation to 
infinite lead length 

As already discussed in Section[21 the limit of an infinitely 
long lead is required in order to obtain the conductance. 
While for the persistent current, we had been restricted 
to rather small ring sizes, the charge stiffness allows us 
to numerically treat rings almost an order of magnitude 
larger. This will enable us to take a closer look at the 
scaling of the charge stiffness with 1/L, even in difficult 
cases like in the presence of transmission resonances. 

As the derivation of the charge stiffness as a function 
of the transmission amplitude in appendix [S] shows, the 
charge stiffness for large rings can be expanded in powers 
of 1/L. In the limit L — > oo, a non- vanishing contribution 
given by (JSJ allows us to determine the conductance. 

The leading corrections lj26|l and Ij29(l for odd and even 
number of particles, respectively, are of order 1/L. Essen- 
tial for the relevance of these corrections is their depen- 
dence on the derivatives with respect to k of the trans- 
mission \t\ and the relative phase shift 5a characteriz- 
ing the scattering region, dda/dk is proportional to the 
Wigner delay time [2E]- At resonances, the two deriva- 
tives may become very large. Then, only rings of circum- 
ference L ^> d|f|/dfc, d5a/dk allow to perform a reliable 
extrapolation to the asymptotic limit. This situation will 
be discussed in Section [3.21 Outside resonances, we found 
that the extrapolation can usually be performed with rings 
about three or four times as large as the scattering region. 

In Fig. fJI we present the deviation of the logarithm of 
the charge stiffness from its asymptotic value as a function 
of the inverse circumference L of the ring. This plot is 
the analogue of Fig. [21 where the scaling of the persistent 
current was depicted, but now the size of the nanosystem 
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is up to a factor of three larger. In all cases shown here, 
we are far away from any resonance. 

The scaling with the ring size has been described by 
different laws in the literature. A parabolic fit was as- 
sumed in the first paper of Favand and Mila JUj while 
a linear fit to the deviations of the logarithm was em- 
ployed in our previous paper |12| . Different polynomial 
scalings were compared by Meden and Schollwock |KS| . In 
the present work, we have used a linear scaling for the 
deviations of ln(D). A second-order fit becomes necessary 
when numerical limitations prevent us from attaining suf- 
ficiently large ring sizes as it has been the case for the 
persistent current (cf. Fig. El - 

For the extrapolation of the charge stiffness in the cases 
presented in Fig. 0] it is sufficient to use the scaling law 

D(U, L S , L) = D^U, Lb) exp ) (10) 

to determine the asymptotic value DoaiU, Lg). The con- 
ductance is then obtained from @ as 

«--•(§£)■ <"> 

This procedure had been used in Ref. J2| to compute the 
influence of the interaction strength on the conductance 
of correlated nanosystems at half filling. The conductance 
of a clean system decreases with the interaction strength 
(see the solid line in Fig. [BJ for even numbers of particles, 
and remains perfect (g = 1) for odd numbers of particles 
independently of the interaction strength. 

3.2 Scaling close to transmission resonances 

The leading correction i|26|) or l|29|) to the charge stiff- 
ness may play an important role close to transmission res- 
onances, where the Wigner delay time and d|t|/dfc are 
large. We illustrate the difficulties in the extrapolation 
procedure present in this case by considering a nanosys- 
tem separated from the auxiliary lead by two tunnel bar- 
riers (cf. Fig.[3J). In order to tune the Fermi energy of the 
ring to a resonance, we introduce an electrostatic potential 
Vo between the tunnel barriers of height Vt> = 1 . A single- 
particle term Vb(ni + riL s ) + Vq J^i^ 1 ni 1S tnus added 
to the Hamiltonian @. The electron-electron interaction 
is present on all Ls sites including the two barrier sites 
but the lead remains non-interacting as usual. We note 
that the additional potential Vo will change the electron 
density in the nanosystem. 

Resonances occur whenever the ground state energies 
of the ring with N + 1 particles and N particles inside 
the double-barrier system are identical. In the absence of 
electron-electron interaction, this implies that the energy 
of the first unoccupied single-particle state of the well be- 
tween the barriers lines up with the Fermi energy of the 
leads. When the degeneracy between ground states with 
different number of particles in the system appears, the 
energetic cost for transporting a particle through the sys- 
tem is zero and the transmission is one. 




Vo 



Fig. 5. Sketch of the site potentials used for the double bar- 
rier system. Electron-electron interaction is present only on the 
grey sites. 




0.0 0.01 0.02 0.03 0.04 0.05 



l/L 

Fig. 6. Scaling towards the asymptotic value of the stiffness 
D for a weakly coupled nanosystem with Ls = 10 and U = 1. 
The circles, squares, and triangles correspond to electrostatic 
potentials Vo = —0.8 (out of resonance), Vo = —1.4 (just to 
the right of a resonance), and Vo = —1.5 (just to the left of a 
resonance), respectively. 



For the reasons discussed in the previous section (see 
also AppendixE}, this case is characterized by a slow con- 
vergence towards the limit Ll — > oo. Large lead lengths 
are then needed because the very rapid changes of the 
transmission as a function of k lead to large correc- 
tions. Another reason consists in the difficulty to main- 
tain the resonance condition for the electron density of 
the nanosystem in the scaling procedure. However, even 
in this unfavorable case, the conductance can be obtained 
by going to larger systems and taking the asymptotic value 
with a greater care than for the non-resonant case. 

Fig. El shows for the example of a double-barrier sys- 
tem how one can extrapolate to the asymptotic value of 
the stiffness in three cases, one favorable and two unfa- 
vorable. The ratio ln(L> / Z?oo) is given as a function of the 
inverse total length of the ring. The circles correspond to 
Vq = —0.8 and U =1, situated in the valley between 
two resonances where the conductance is small. In this 
case, the extrapolation is straightforward and the slope 
is very small. The other two cases are different. Taking 
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Vq = — 1.4 and [7 = 1 (depicted by squares), we are just 
to the right of a resonance. The corrections to the scal- 
ing formula i|10|) are very large for small ring sizes, and 
a naive extrapolation from there can give wrong values 
(even g > 1) for the conductance. In order to test that the 
asymptotic value for D is approached, one calculates the 
parameters C and of the scaling formula l|10JI for two 
different values of L and one continues to increase L until 
the slope C and the asymptotic stiffness Doo converge to 
constant values. In the case Vq = —1.5 shown by triangles 
in Fig. we have first determined C and assuming 
the scaling law (|10|l for L = 30 and L — 40. Because the 
procedure gives different results when we take L = 40 and 
L = 50, we were forced to increase L. Since the values 
for C and obtained with L = 120, 130 and 140 do 
not vary, we assume that one has reached the asymptotic 
regime. This procedure can require large values of the to- 
tal length L of the ring, which arc difficult to reach for 
large filling factors v. Using a fit with more parameters 
can be an option when the convergence is slow, but the 
extrapolation must be done very carefully. The behavior 
of the stiffness D as a function of the length L in this last 
example is quite complicated because the density in the 
lead cannot be kept perfectly uniform and therefore the 
resonances move as a function of the increasing size of the 
ring. This extreme case illustrates the potential difficul- 
ties which must be solved in order to get reliable values 
for g in the vicinity of transmission resonances from this 
method. For Vb = —1.5 (just to the left of a resonance), 
the slope has changed sign and we still need to go to big 
ring sizes for a reliable extrapolation. 

In Fig. CK we depict the results of the conductance, 
evaluated using the previous extrapolations for the two- 
barrier system. We compare the results for [7 = and 
[7=1. The values for U = have been obtained in the 
same way as the values for U = 1, using DMRG and the 
extrapolation. They are found to agree with results from a 
non-interacting Green function calculation. The fact that 
we do obtain perfect conductance (g = 1) at resonances 
supports our claim that the asymptotic procedure is capa- 
ble of yielding the correct transport properties. In Fig.Q:> 
we show the slope C(U,Ls) of the scaling law (jTUj) . As 
one can see, the resonance structure is clearly reflected by 
the slope of the scaling curves. The jumps in the slope 
coincide with the values for which the dimensionless con- 
ductance approaches its maximum value of one. The slope 
is closely related to the behavior of d|i|/d/c. As expected, 
the interaction U changes the position of the peaks and 
their widths. 



4 Conductance outside half filling 

As stated in the introduction, most of the applications 
of the embedding method have so far been restricted to 
half filling. In the previous section, we have maintained 
half filling in the average over the composed ring, but the 
filling of the correlated system itself depended on the po- 
tential Vo between the barriers. As an even more general 
situation, we now consider the case where the filling in 




-5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 



V 

Fig. 7. (a) Conductance g and (b) slope C(U,Ls) of the 
scaling law IIUI are shown as a function of the electrostatic 
well potential Vb for the configuration of Fig. |S| Results for 
interaction strength U — 1 are indicated by full symbols and a 
solid line while the open symbols and the dotted line represent 
results for the non-interacting case (U = 0). 



the composed ring has an arbitrary value v. The half-filled 
systems exhibit particle-hole symmetry, and therefore the 
compensating potential V+ = 1/2 required to yield charge 
neutrality inside the nanosystem is known a priori. If we 
want to ensure a given constant filling v for the nanosys- 
tem and the lead even when the lead length is changed, V+ 
becomes a function of the interaction strength and the ring 
size. In this section, we extend the method to nanosystems 
outside half filling which are well coupled to the lead. By 
choosing the appropriate particle number, the same fill- 
ing is imposed in the auxiliary lead in order to obtain the 
transmission coefficient \t(E-p,U)\ 2 at the corresponding 
Fermi energy and to ensure a better convergence towards 
the limit of infinite lead length. 

In order to determine V+ for an arbitrary filling v, 
we begin with an initial guess for V+ and calculate nu- 
merically the corresponding number of particles contained 
inside the nanosystem. Then, we adjust V+ performing 
an iterative solution of the problem using the Newton- 
Raphston method. In principle, V+ depends on U and L^. 
For example, for Ls = 8, v — 3/8 and U = 3, V+ varies 
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Fig. 8. Conductance as a function of U for different filling 
factors v of a correlated system of size Ls = 8. 

from 0.1924 to 0.1939 as L L is doubled from 24 to 48. At 
a fixed interaction strength, the dependence of V+ on 
becomes negligible beyond a certain Ll, and can then be 
ignored. Therefore, the iterative procedure has only to be 
performed until a limiting value for V+ has been attained. 
Then this value can be kept for larger ring sizes from which 
C and Dec are determined, using the same scaling law as 
at half filling. 

In Fig. the conductance of a nanosystem of length 
Ls — 8 perfectly coupled to the lead is given as a function 
of the interaction strength U at different filling factors v. 
Since the filling is kept uniform everywhere in the ring, the 
curves characterize g(E-p,U) at the corresponding Fermi 
energy Ep. At v = 1/8 (short dashed line), in average 
only one particle is left in the nanosystem. In the absence 
of other particles to interact with, the dimensionless con- 
ductance therefore equals one, independently of the inter- 
action strength. For larger filling factors, the conductance 
g decreases with increasing interaction strength U and 
this decay becomes more pronounced as the filling factor 
is increased. The rather sharp drop of the conductance 
occurring at half filling around U = 2 is a precursor of 
the Mott transition expected in the thermodynamic limit. 
The conductance above half filling can be obtained from 
g(v) = g(l — v) as a consequence of particle-hole sym- 
metry. The influence of the interaction strength on the 
conductance is thus the strongest at v — 1/2 as expected. 



5 Conductance for disordered nanosystems 

Having demonstrated that the conductance of a correlated 
nanosystem can be obtained from the charge stiffness af- 
ter embedding it into a large noninteracting ring, we now 
apply this method to the problem of interacting electrons 
in disordered systems. The effect of repulsive interactions 
in a disordered system is a controversial issue - It is 
often believed that interactions impede transport. This 
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Fig. 9. Scaling towards the asymptotic value of the stiffness 
D for disordered samples (W = 5). For the same disorder re- 
alization, two values of the interaction are shown for half fill- 
ing. Open triangles represent (7 = 4 and open circles the case 
[7 = 1. The negative slope in the former case corresponds to a 
charge reorganization and the conductance g = 0.34 for U = 4 
is greater as compared to g = 0.018 for (7 = 1. Results for 
the same disorder configuration with (7 = 4 and v = 3/8 are 
displayed with filled triangles. The conductance in this case 
is g = 0.0011, demonstrating that the charge reorganization 
depends on the filling. 

belief comes from perturbative arguments showing that 
interactions reduce the density of states at the Fermi level 
of a disordered metal jSH] and open a gap for a strongly 
disordered insulator . On the other hand, in the strong 
disorder limit zero temperature transport can be enhanced 
by an interaction- induced delocalization of the many-body 
ground state. This was demonstrated for the special case 
of half filling in Ref. fT^]. In the following, we will study 
the role of the filling factor in the delocalization process. 

We include the disorder potential into the Hamiltonian 
(TJJ by adding a term 

Ls 

H dis = Wj2v ini , (12) 

i=l 

where W denotes the disorder strength, and the u$ are in- 
dependent random variables, equally distributed within 
the interval [—1/2,1/2]. The disorder potential is only 
present within the nanosystem of length Lq. 

We start by verifying that the scaling towards infinitely 
large rings also works in the presence of disorder. Fig. [5] 
depicts the dependence of the logarithm of the charge stiff- 
ness D on the ring size L for a sample with W = 5 for 
interaction strengths [7=1 (circles) and 4 (triangles). The 
disorder realization is the same in both cases. The open 
symbols refer to v = 1/2 while the full symbols correspond 
to v = 3/8. In all cases the scaling works well, and thus 
reliable values for the conductance can be extracted. 

The analysis of individual samples helps us to under- 
stand the physical mechanisms involved when disorder 
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Fig. 10. Logarithmic ensemble average of the conductance as 
a function of the nearest neighbor repulsion U for a disordered 
nanosystem of length L$ = 8. The open symbols correspond 
to half-filled nanosystems, and the filled symbols to a filling 
factor v — 3/8. The triangles correspond to W = 1, the circles 
to W — 5 and the diamonds to W = 9. For W = 0, the dashed 
line corresponds to v = 3/8 and the dotted line to v = 1/2, 
respectively. 

and interactions are both relevant [H2HSB] - Studying the 
evolution of the ground state energy or the electron den- 
sity as a function of U, we can detect charge reorganiza- 
tions at critical values of the interaction strength. For the 
sample shown in Fig.UHwe have, at half filling, a charge re- 
organization in the ground state structure around (7 = 4. 
Charge reorganizations appear when a ground state con- 
figuration which is well adapted to the non-interacting 
case, where the fermions are located in the minima of the 
disorder potential, changes towards a Wigner-like crys- 
talline structure which is energetically favorable at strong 
repulsive interaction. This resonant situation increases the 
conductance at the particular (sample dependent) cross- 
over value of the interaction. In other samples the charge 
reorganizations can occur at different values of the inter- 
action or can even be absent, depending on the disorder 
realization. Reducing the filling makes the charge reor- 
ganizations less likely. For the charge reorganizations of 
the disordered case, we typically obtain a negative slope 
for the asymptotic scaling of D, similar to the case of 
clean systems with odd number of particles [TjJ. In both 
cases, a degeneracy of different charge configurations in 
the nanosystem occurs. 

In Fig. ^3 the ensemble average of the logarithm of g 
is given as a function of U for disorder strengths W — 1 
(triangles), W — 5 (circles), and W — 9 (diamonds) and 
filling factors v — 1/2 (open symbols) and v = 3/8 (full 
symbols). One can see from the increase of the average 
conductance at weak interaction in the strongly disordered 
case, W = 9, that the nearest neighbor interaction has 
stronger derealization effects around half filling. The re- 
sults for the clean case show Mott insulator like behavior 
at half filling (dotted line) . The decay of the typical value 



of g as a function of U is faster for the half-filled case 
than for v = 3/8. When we introduce a random potential 
in the nanosystem, the reduction of the typical conduc- 
tance due to localization effects is more important outside 
half filling. 

The larger the density, the better is the screening of the 
random potentials. In the case of weak disorder, W = 1, 
this gives rise to a crossing of the curves with v = 1/2 
and v = 3/8 as U increases. For stronger disorder, this 
crossing occurs at larger values of U (U = 4.5 for W — 5) 
and for very strong disorder (W = 9) the crossing cannot 
be observed in the figure. 

We can also see in Fig. ED that in the strong disor- 
der case (here W = 9), nearest neighbor interactions can 
favor transport. This enhancement of the typical elastic 
transmission, and hence of the zero temperature conduc- 
tance, is maximal around U = 0.5 and, though mainly 
characteristic for half filling, it persists outside v = 1/2. 

The charge reorganization induced by repulsive in- 
teractions in strongly disordered systems and its associ- 
ated derealization effect was first observed in the persis- 
tent current of nanosystems 32 , 33 forming a ring (with- 
out the auxiliary lead introduced within the embedding 
approach). As our results demonstrate, the same effects 
can be found in the conductance g. Considering a given 
nanosystem, one observes a similar resonance structure 
[T^] as for the persistent current |32U33j . although the in- 
dividual peaks are wider for the conductance than for the 
persistent current. 



6 Summary 

The residual conductance of a correlated nanosystem can 
be obtained from the charge stiffness or from the persistent 
current of a ring composed of the system and an auxiliary 
non-interacting lead. Using DMRG for spinless fermions, 
we have numerically studied basic properties of this em- 
bedding approach. In particular, we have demonstrated 
that the flux dependence of the persistent current for an 
interacting system and a non-interacting lead agrees with 
the flux dependence of a non-interacting ring with a scat- 
terer, in the limit of infinite lead length. This allows to 
extract the interaction dependent transmission coefficient 
of the interacting system, and hence its residual conduc- 
tance. 

A detailed analysis of the finite-size corrections has 
been performed for the charge stiffness. The main features 
of these corrections can be understood from the analysis 
of the non-interacting case. Away from transmission res- 
onances, we obtain a very good scaling behavior already 
for not too large lead lengths, and the conductance of the 
correlated nanostructure can be readily obtained. Close 
to resonances the asymptotic limit of large lead lengths 
is problematic and only by considering very long leads 
we obtain the correct asymptotic behavior. Even in these 
special cases, the results for the conductance agree with 
our expectation for the resonant tunneling behavior in a 
double barrier structure. 
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It is straightforward only at half filling to keep the elec- 
tron density in the correlated system fixed when chang- 
ing the ring size. We have demonstrated that an exten- 
sion of the method to arbitrary filling factors is feasible 
provided the compensating potential is adjusted appropri- 
ately. For clean samples, it was observed that the decrease 
of the conductance with increasing interaction strength is 
strongest at half filling and becomes weaker as the filling 
factor changes towards smaller or larger values. 

Another extension consists in the introduction of dis- 
order in the correlated system. Charge reorganizations of 
the ground state appear at sample-dependent values of 
the interaction strength, affecting the long lead scaling 
and the asymptotic values. In the ensemble averages, we 
obtain for weak disorder a decreasing conductance as a 
function of the interaction strength. However, for strong 
disorder we have shown that a nearest neighbor repulsion 
can enhance the average of the logarithm of the conduc- 
tance for spinless fermions in one-dimensional samples. 
This enhancement persists outside half filling, although it 
becomes weaker. 

So far, the approach is still limited to spinless fermions 
and single-channel leads, although the system itself can 
be arbitrary. Nevertheless, the method is well suited to 
study the role of the contacts between the nanosystem 
and the leads. Furthermore, interesting phenomena like 
even-odd oscillations of the conductance with the num- 
ber of fermions were found with this approach |12II34| . In 
the absence of spin-flip scattering, the generalization to 
electrons with spin is straightforward. Indeed, first cal- 
culations for the Hubbard model have already been per- 
formed |12j . These and further issues will be explored in 
more detail in future work. 

RAM wishes to thank J. Segala for reminding him of some 
properties of the Chebyshev polynomials. We gratefully ac- 
knowledge financial support from the European Union through 
the RTN program (Contract No. HPRN-CT-2000-00144). PS 
was supported by the Center for Functional Nanostructures of 
the Deutsche Forschungsgemeinschaft within project B2. 

A Flux dependence of the ground state 
energy for large rings with a small 
non-interacting scattering region 

In this appendix, we discuss the flux dependence of the 
ground state energy for a ring containing a non-interacting 
scatterer. The scattering region of length Lg is connected 
to a disorder-free lead of length L L • This arrangement is 
closed to a ring of total length L = L§ + L L , as shown 
in Fig. ^ We present a systematic expansion in powers 
of 1/L starting from the limit of infinite lead length for 
the flux-dependent part of the ground state energy. This 
leads to analytic expressions for the asymptotic values of 
the persistent current and the charge stiffness, as in Ref. 
|25j . We extend this theory by calculating the first finite- 
size corrections to the flux-dependent part of the energy 
and the charge stiffness. These corrections are important 



to understand the way in which the asymptotic values are 
approached when we extrapolate to infinite ring size. 

The one-particle eigenenergies of the ring are given by 
the quantization condition 

det (I - Ml Ms) = , (13) 

where Ms and Ml are the transfer matrices of the system 
and the lead, respectively. In the presence of time-reversal 
symmetry, the transfer matrix of a one-dimensional scat- 
terer can be expressed in terms of three independent an- 
gles a, 9 and tp: 

/ • x (14) 

1 / e m /sin# — i cot 9 + cos tp \ 

shop \ i cot 9 4- cos tp e~ iQ /sin# J ' 

where the two components correspond to right and left 
moving particles while r and t are the reflection and trans- 
mission amplitudes, respectively. The angle a is the phase- 
shift associated with the scattering region. Whenever the 
right- left symmetry is respected, we can set tp = x/2, and 
the expression of Ms simplifies considerably. However, this 
symmetry requirement is not satisfied for disordered sam- 
ples. In the general case the transmission amplitude is 
given by t — e 1Q sin 9 sin tp. 

The transfer matrix of a lead of length Ll for a state 
with wave number k > reads 

Here, we have made use of the fact that the flux can be 
transformed into a boundary condition which may be pre- 
scribed in the lead. 

Inserting the transfer matrices (|14J) and ljl5|l . the eigen- 
value condition (|13l) yields 

cos(<£) = yj— cos (kL + Sa(k)) . (16) 
\t(k)\ 

Here, we have introduced the phase shift 5a — a — kLg 
of the scattering region relative to a perfect lead of the 
same length L$- The solution of (TTHjl yields the quantized 
momenta k of the energy eigenstates in the lead. 

Since both, t and 5a are functions of k, it is in general 
impossible to obtain an analytic solution of (|16fl . However, 
progress can be made in the asymptotic limit of large L, 
which was worked out by Gogolin and Prokof'ev |2S] in 
their study of the persistent current. We extend their ap- 
proach to calculate the first finite-size corrections of the 
charge stiffness. Furthermore, a generalization to arbitrary 
dispersion relation in the lead allows us to discuss contin- 
uum and tight-binding models at the same time. 

The eigenvalue condition (|16f) can be rewritten as 

k = k° + jf ± (k,ci>) . (17) 

Here, k^ — 2irn/L with n > denotes the eigenvalues in 
the case of perfect transmission with |i| = 1 and 5a = 0. 
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Following the notation of Ref. |25| . we have furthermore energy then reads 
introduced 



f±(k,$) = ±Arccos(|t(fc)| cos<P) - 8a(k) 



(18) 



By Arccos, we denote the principal branch of the inverse 
cosine function that takes values in the interval [0, it]. In 
order to ensure a positive value for fc, f^(k,<P) should not 
be used for the case n = 0. The splitting of the solu- 
tions of ifTTjl corresponding to "+" and "-" cannot exceed 
the spacing 2ir/L between the k^, provided that Sa(k) is 
smooth on this scale. This is the case in the limit L — > oo 
and ensures that the order of the solutions with respect 
to energy is given by n. 

Iterating l|17l) and expanding f± for large systems, we 
obtain the expansion 



df±(k,< 



1 d 
2lJdk 



9k / k=k o 

df±(k,<py 



Ok 



o 



k=k<± 



(19) 

for the solutions of (|17|) in powers of l/L. Such an ex- 
pansion is problematic in the vicinity of resonances, when 
dSa/dk and d|t|/dfc are very large. Then, the expansion 
is valid only for sufficiently large L. 

We now calculate the ground state energy of the sys- 
tem as a function of the flux to order l/L 2 . The dispersion 
relation in the perfect lead will be denoted by e(k). Using 
H19(l . we start by expanding the one-particle energies in 
powers of l/L and obtain 



E^i^^eik^ + J^Kk+J + eik-)} 



= + (|p [Arccos(|t| cos<£) - 6a] 2 



k=Q 



5a{k%) 



k=kO 



+ [A 1 'ccos 2 (| i |cos^) + ^])^ o 



1 d 2 
ZL^dk 2 



de 
Ok- 



[3foArccos 2 (|t| cos<P) + 8a 3 



k=k° 



(21) 

The sum runs up to rip = (N — l)/2. We have assumed 
(de/dk)k=o — and kept all terms which can give rise to 
contributions up to order l/L 2 . The first term in the sum 
is the ground state energy in the absence of scattering. For 
finite filling, i.e. for N of order L, it is proportional to L 
while the second term representing the energy change due 
to the scattering potential is of order 1. The third and 
fourth terms are the leading flux-dependent corrections. 
Since we are interested in the persistent current and the 
charge stiffness, these are the only terms in the sum which 
need to be considered further. Converting the sums over n 
into integrals, these flux-dependent contributions can be 
expressed as 



+ 



L \dk 

1 d (de 
'Ztfdk \dk 

1 d 2 /de 



6L 3 dk 2 V dk 



f±(k,$) 



k=k° 



k=k» 



k=k° 



o 



L 4 



(20) 



For an odd number N of spinless electrons in the ring, 
all occupied states n come in pairs ([n,-] and [«,+]), except 
for the one corresponding to n = 0. The total ground state 



kp+TT J L 



2^L I dfc|-(|Arccos 2 (|t|cos<P) 

tt/L 

= ^pArccos 2 (|i(fc F )| cos(0)) 



dk 2 



Arccos 2 (|t| cos^) 



k=0 
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and 



(|21[1 and obtain 



1 f d 2 ( de \ 
l^L 2 J M—^-6aA rC co^\t\co S <P)) 

-<5aArccos 2 (|i| cos <P) 
- ( ■^4fcArccos 2 (|i| cos^)^) 1 




2ttX 2 1 dk\dk 

}2, 



"'Hi- 



(23) 



respectively. Here, fcp = 2im-p/L is the Fermi wave num- 
ber and vf = (de/hdk)k=k F is the Fermi velocity. Taking 
the derivative of the leading flux-dependent term of the 
ground state energy 



E 



odd(i; 




2^L 



Arccos 2 (|i(fc F )|cos(<?)) (24) 



with respect to the flux cj>, one obtains the asymptotic form 
of the persistent current given in Q for an odd number of 
particles. The leading order of the charge stiffness of ||HJ| 
is obtained as 



D 



(i) - 



L 



(£° dd(1) (0) 



- e: 



,odd(i; 



(7T) 



Arccos(|i(fc F )|) 



(25) 



As we will show below, this last result is independent of 
the parity of the number of particles. 

The first finite size-correction to these asymptotic val- 

(2) 

ues follows from the second-order contribution Eq to the 
total energy. Using the terms of order 1/L 2 from l|22ll and 
<|23|) . and taking into account the contribution from the 
particle in the state [0,+] in the second line of i|2 ljl. we 
obtain the correction to the charge stiffness for an odd 
number N of particles 



£) odd(2) 
1 

~~ ~ 2L 



L 



e: 



odd(2) 




(0) 



E° dd(2) (7T) 



hvp(Sa — 7r) 



dSa 
"dk 



d\t\ 



dk 



— — Arccos(|t| 



k—kp 



dk 2 



Arccos( 

2 v 



k=0 



(26) 

The last term vanishes if we assume that \t(k = 0) j = 0. 

In order to treat also the case of an even number of 
particles, we subtract the contribution of the particle in 
the one-body state [np, +] from the total energy of Eq. 



E° dd (<P)-e(k F )~ 



1 f de 



L\dk 



1 9 f de *2n ^ 



k—kfT 



k—kp 



(27) 

With these additional terms, one obtains the leading 
flux-dependent term of the ground state energy for an even 
number of particles as 



E 



cvcn(l) 



2ttL 



Arccos 2 (|i(fc F )| cos(# - tt)) . (28) 



The derivative with respect to <p leads to the asymptotic 
form of the persistent current of Eq. J5J for an even num- 
ber of particles, and the result for the leading contribution 
to the charge stiffness agrees with (I25|) . 

For the first finite-size correction to the stiffness we 
obtain 



rjeven(2) Zl ( . 
1 

~ ~ 2L 



H2) 



1(2) 



(TT) 



d 2 e 



dSa 

— Sa + hvF — 



(f -Arccos(|i|) 



(29) 



hvpSa 



d\t\ 



1 



dk ^/T^W 



k—kp 



which differs from the case of an odd number of particles 

^Frorn Eqs. J23) and (j^U one can see that the 1/L 
scaling for approaching the asymptotic values of the stiff- 
ness is problematic close to resonances, where d5a/dk 
and d|i|/dfc are large, and \t\ approaches 1. Assuming an 
isolated Breit-Wigner resonance the Wigner time is 
proportional to g and the corrections L>' 2 ' are essentially 
given by the half width of the resonance. Outside reso- 
nances where 6a -C 1 and for small \t\, the leading correc- 
tion to the stiffness can be approximated by 



D 



(2) 



1 . .d6a 
2L hVFltl ^k 



Therefore, one obtains for this case 



In 



D 



In 



1 d6a 
L~dk 



(30) 



(31) 



and the Wigner time gives the slope of the scaling curve. 
The above arguments are valid in the non-interacting case. 
However, the intuition developed in this case is also useful 
to interpret our numerical results for the interacting case. 



B Conductance of a NSN region from 
persistent current 

In this appendix we treat the case of a superconductor 
between two metallic leads. This is a striking example of 
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a correlated system exhibiting non-Fermi liquid behavior. 
It will be demonstrated that the correct result for the 
conductance can be obtained from the persistent current 
by means of 



B.l Double Andreev scattering 

The Andreev scattering at a NS junction, i.e. the interface 
between a normal metal and a superconductor, is a well- 
known phenomenon. In an Andreev scattering process, an 
electron coming from the normal metal is reflected as a 
hole while a Cooper pair moves on in the superconductor. 
The linear conductance of the interface between the nor- 
mal metal and the superconductor in the one-channel case 
is given by 

4e 2 T 

(32) 



G = 



h 2-T 



where T is the transmission probability in the normal 
metal [3l)ll3*7] . For the normal lead, T = 1 and one gets 
that the resistance of a single normal-superconductor in- 
terface is the half of the resistance without interface. 

In the following, we will consider a NSN double junc- 
tion consisting of a clean superconducting layer of thick- 
ness Lg connected to normal-metal electrodes by perfect 
interfaces. It is assumed that the superconducting gap 
A{x) jumps at the interface from zero in the normal metal 
to its full value A inside the superconductor 



A{x) = AO(x)0(L s - x) , 



(33) 



where 0(x) is the step function. This approximation is 
common in the treatment of mesoscopic superconductors 
|38|. Blonder, Tinkham and Klapwijk calculated the 
conductance by solving the Bogoliubov-de Gennes equa- 
tion with this rigid-boundary condition and found for 
T = 1 the linear conductance G — 2e 2 /h. This result 
can be understood by taking two Andreev interfaces with 
conductance (|32|l in series. 

When we close the two normal metal leads of the NSN 
junction to a ring, we recover the geometry of the embed- 
ding method where the correlated system is formed by the 
superconductor. It is therefore interesting to see how one 
can recover the linear conductance from this approach. 



B.2 Persistent current and conductance of a 
superconductor 

For a one-channel ring consisting of a normal conduct- 
ing region of length Ln and a superconducting region of 
length L$, the solution of the Bogoliubov-de Gennes equa- 
tion for a boundary condition analogous to 1331) yields the 
persistent current |4()U41| 



m = - 



ev-p 



n L N + £ tanh(L s /£ ) 



ET m {X) . 
sin(TO$>) 



Here, £o = hvp/A is the superconducting coherence length 
and T m (X) denotes a Chebyshev polynomial in the vari- 
able 

cos(fc F L) 



In the limit £o - 
ring of length L 



cosh(L s /£o) 
oo, one obtains a normal conducting 
Ln + Lg with the persistent current 



jo ( <£) = 1 y I r sin ( m (,|> _ kpL) \ 

m— 1 



(36) 



sin (m(<£ + k F L))] 



Apart from a factor of two accounting for the spin, this 
expression reduces to Q or |J5} for \t\ = 1 depending on 
the parity of the number of particles per spin. We note, 
however, that the expression l|34|) for the NS ring can, in 
general, not be expressed in the form or J5J. 

According to @, the dimcnsionlcss conductance g can 
be obtained from the persistent current at flux <P = tt/2 



J(tt/2) 



7T L N +£ tanh(L s /£o) 



— |sin m ^— — Arccos(AT) 



+ sin m ( — h Arccos(AT)^ j . 

2 (37) 
By means of the Fourier representation of a sawtooth func- 
tion, one finds that the absolute value of the persistent 
current becomes 

efF 



J(tt/2) 



L N + f tanh(L s /£o) ' 



(38) 



In view of this result, the superconducting region can be 
thought of as a normal-conducting metal of an approx- 
imate effective length given by the minimum of Lg and 
&• 

It is now straightforward to determine from l|38(l the 
dimensionless conductance 



lim 



lim 



J°(7r/2) 

£n + £s 

£n + £o tanh(Ls/£o) 



(39) 



(34) 



Here, the persistent current of the normal ring can again 
be thought of as being obtained from l|38|) in the limit 
£o — * oo. We thus recover the correct result g = 1 for 
the dimensionless conductance. The leading corrections 
depend on the ratio [Lg — Co tanh(Ls/Co)]/-^N between 
the relative length of the superconductor, i.e. the differ- 
ence between the real length of the superconducting re- 
gion and its effective length, and the length of the normal 
region. Even though here the transmission amplitude re- 
mains equal to one in the presence of correlations, this ex- 
ample gives another demonstratation that the embedding 
methods works, even for having the conductance through 
a system which is very far to exhibit a Fermi liquid be- 
havior. 
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